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The nonasymptotic critical behavior of sound attenuation coefBcient has been studied in an 
elasticaUy isotropic Ising system above the critical point on the basis of a complete stochastic 
model including both spin-energy and lattice-energy modes linearly coupled to the longitudinal 
sound mode. The effect of spin-lattice relaxation on the ultrasonic attenuation is investigated. 
The crossover between weak-singularity behavior and strong-singularity behavior ig 

studied as dependent on the values of ultrasonic frequency, reduced temperature, relaxation times 
etc.. A new high-frequency regime with a singularity of the type t~^'^^°' is discovered in the magnetic 
systems. This new regime corresponds to an adiabatic sound propagation and is very similar to the 
ones in binary mixture and liquid helium. The scaling functions are given in all regimes to the 
first order in e. A new frequency-dependent specific-heat being the harmonic average of the bare 
lattice and critical spin specific-heats is introduced. It was shown that such specific-heat describes 
the process of equilibrization between spin and lattice subsystems and includes the most important 
features of critical sound attenuation. In some regions of coupling constants the acoustic self-energy 
can be very well approximated solely by this quantity. 
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I. INTRODUCTION 

A strong anomaly of sound attenuation measured in Ising like magnetic metals [^j^ is well described by the theories 
[H^illljll ill which the sound wave is assumed to couple to two critical spin fluctuations. The sound attenuation 
coefficient is found then as the imaginary part of the four-spin response function. On the contrary to metals, in 
magnetic insulators we observe a weak anomaly in sound attenuation [|l],P, which has been phenomenologically 
explained by postulating a linear coupling of the sound mode to the spin-energy with the spin-lattice (or more 
correctly spin-energy-lattice-energy) relaxation time playing the essential role. Assuming only the specific-heat like 
singularity of the spin-lattice relaxation time and simple separability of contributions coming from both types of 
coupling, Kawasaki 0| was able to obtain weak divergence of the sound attenuation coefficient proportional to the 
square of the specific-heat for the magnets where the first coupling is negligible. However, as it was shown by 
the renormalization-group analysis of the model of coupled spin and energy fields |^ the first assumption is not 
correct sufficiently close to the critical point where the energy relaxes with the same characteristic exponent as the 
spin fluctuations, ze = because the singular part of the energy response function is dominated by the four-spin 
response function. Also the second assumption concerning the separability fails near Tc as was recently shown 
The objective of this paper is to use a complete stochastic model which includes interactions between the order 
parameter (spin), the acoustic phonon, spin-energy and lattice-energy modes in order to find the asymptotic as well 
as nonasymptotic behavior of the sound attenuation coefflcient. We shall restrict our discussion to the disordered 
phase. We consider an elastically isotropic Ising model, where above the critical temperature a phonon (longitudinal) 
mode Q is coupled to a scalar order parameter S as well as to the fluctuations of spin- and lattice- energy densities 
es and e^. Introduction of the second energy density associated with the lattice (or conduction electrons in some 
metals) allows us to consider much richer dynamics than that of the two limiting cases: pure relaxation to a reservoir 
with inflnite thermal conductivity or pure diffusion of energy. Our model permits also a study of the intermediate 
dynamics between these limiting cases. We have obtained an general expression for the sound attenuation coefficient 
properly describing the ultrasound attenuation in very broad temperature and frequency range. The purpose of this 
paper is to demonstrate that the strong sound attenuation anomaly (sometimes called Murata-Iro-Schwabl behavior) 
as well as the weak anomaly (Kawasaki behavior) can be obtained within our stochastic model depending on the 
relative size of the reduced temperature, frequency and other parameters of the model. The nonuniversal amplitudes 
for both types of behavior are modified considerably by the presence of the second energy field in comparison with 
one-energy-fleld model. In the high-frequency regime we obtain a new kind of critical behavior, which we shall call 
the "adiabatic" behavior, analogous to the critical attenuation behavior in ^He and in binary mixtures. Contrary 
to the Murata-Iro-Schwabl and Kawasaki behavior the amplitude for the "adiabatic" limit is determined completely 
only by one coupling constant, the one describing the interaction of the sound mode with two fluctuations of the 
order parameter. The "adiabatic" limit does not change by the inclusion of the field ez,. We shall also show that in 
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some regions of parameters the acoustic self-energy is simply proportional to a frequency-dependent specific-heat C_ 
which is the harmonic average of the bare lattice specific-heat and the specific- heat of an idealized spin system (model 
A in terminology of C_ has a simple interpretation : it equals to the ratio of the heat transferred from one 

sub-system to the other, to the induced temperature difference between subsystems. The important point is that such 
specific-heat shows three main types of singularities characterizing the sound attenuation coefficient as the critical 
temperature is approached for low or high frequencies. 

The paper is organized as follows. In Sec. II we introduce the model and then perform some dynamical decoupling 
transformations. In Sec. Ill we obtain a general expression for the acoustic self-energy taking into account the re- 
ducibility of the latter with respect to energy propagators. In Sec. IV many different regimes for sound-attenuation 
are discussed as a function of frequency, temperature and relaxation times. The scaling functions are determined 
within the one-loop approximation. We then show that for some regions of parameters the ultrasonic attenuation can 
be written as the imaginary part of C-{uj). The relevance of this specific- heat to the description of nonasymptotic 
sound attenuation is shown. In Scc.V we smnmarize the results. 



II. MODEL AND BASIC FORMULAS 



A. Equations of motion 

The entropy of the system ||l^,|ll| can be written as 



-Q/3 



2 2 

+ 2gY,e^^S'' + 2fesS'' + 2w{es + aeL)Y,e^^ + ^ + ^\ , (2.1) 

a a ) 

where eap{x) are components of the strain tensor and the symbols Cap stand for the bare elastic constants and Cs 
and Cl are the spin and the lattice specific- heat, respectively. The unitary mass density and ksTc — 1 have been 
assumed. The first three terms in the total Hamiltonian- describing the static behavior of the system, make the 
Ginzburg-Landau functional for the spin variable. The elastic energy is given by the 4th and 5th term in Eq. ( ^.l| ). 
Here, we have made use of the relation Cn — C12 = 2C44 applicable to the isotropic cubic systems. The last two 
terms results from lowest order expansion of the entropy functional with respect to energy fields. The other terms in 
the Hamiltonian describe interactions. The constant g denotes the coupling of the (longitudinal) sound mode to two 
spin fluctuations. The coupling of the sound to the energy fields is characterized by the constants w and aw where 
the parameter a is the ratio of the couplings of the phonon mode to the lattice-energy ei, and spin-energy density eg, 
respectively, / is the coupling constant between the order parameter and the spin energy generating the divergence 
of the specific-heat. 

After introducing the normal mode expansion of the strain tensor, the dynamics of the system can be described by 
the coupled Langevin equations |^,|ll| 

f) H 

5k = -r-— -fCk, (2.2a) 
00-k 



^k = -77^ ^fc'Qk + 77k, (2.2b) 



5Q- 



ek = -(7s + Asfc2)j^+7n-+¥'k, (2.2c) 



SU 5H 

ek = -(7L + Asfc')^+7^+V'k, (2.2d) 

where Qk is the longitudinal phonon normal coordinate and Ck, '7k, <^k and "(/ik are the Fourier components of Gaussian 
white noises with variances related to the bare damping terms F, Ok^ ^ (7s + As'fc^) and (7l + AlA:^) through the Einstein 
relations. 
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The first two of these equations have been commonly used in investigation of critical sound propagation ,|^,|l^,|l3[ . 
The other two describe the energy flows between the sub-systems. In the absence of the nonlinear terms and noises 
they transform into the equations describing the decay of total energy e = es + and equalization of temperatures 
of both sub-systems 

e = -^e-AAeC_A (2.3) 



A = -^e-AC_A (2.4) 

where A = 6Tl — 6Ts = sl/Cl — es/Cs is the temperature difference between both sub-systems, A-)_ = "fs + Jl — 
27 + X+P, AAe = - , A = + 2i+^ + _|2_, X+ = Xs + Xl,C+ = Cs + Cl, and 

C_ = (1/C's + I/Cl)"^ . If 7s = 7l = 7, then the total energy is conserved (for k 0) and ^ = is the rate of 

the thermal conduction with As and Xl being the thermal conductivity for spin and lattice sub-systems, respectively. 
For this case, the temperature difference A decays at the rate AC_ w 7/C'_. Therefore, AC_ can be interpreted as the 
frequency of the spin-lattice relaxation. It reduces to 7/Cs for the lattice of infinite specific heat. For non-conserved 
total energy, as e.g. for 7l > 7s = 7 (we assume that the spin-energy can change only through interactions with 
lattice degrees of freedom i.e. 7s = 7 ), 

A+ A+C^ 

e = -^e ^ ' 



A = -^^e-AC7_A (2.6) 

with A-)_ ~ 7i — 7. The total energy, in the absence of the temperature difference A, decays at the frequency 

(7L — 'y)/C+ and the temperature difference relaxes (in the absence of total energy fluctuations) at the rate AC- = 

1 A+Cg 

C- ^ C+Cl ■ 

In a typical case is the sum of energies of all phonon branches except for the longitudinal acoustic phonons with 
wave- vectors restricted to the sphere < |k| ^ 1. These phonons are taken explicitly into account in the model as the 
quantity of our primary interest in this paper is the phonon response function. In some metals, where the lattice can 
be treated fully adiabatically (Cl — > 00), cl may be understood as the energy density of conduction electrons. In that 
case, the total energy in the system composed of the localized electrons and the conduction electrons sub-systems is no 
longer conserved so we have admitted the possibility of nonconserved total energy cl + es {il 7^ 7) in the equations 
of motion. 



B. Lagrangian 



Instead of using Eq. ( |2.2[ ), it is convenient to represent the model in terms of the equivalent functional form [|14 15 
with a Lagrangian given by 



,S- 



k.- 



Ok (5k,wQ-k,-L 
dH 



rs(fc)gLg^k.-. 



(-tj^ -|-i6'fc^u;)(5_k-^ + 



aQk, 

dH 



*t^e^k,-c. 



TL{k)e 



dH 



27ek.^e^k.-. 



dH 



9eL 



(2.7) 



where S,Q,es and cl are auxiliary "response" fields and Ti{k) = 7.; + A^fc^ for i = S,L. The Lagrangian does not 
contain the Jacobian because we have assumed the Heaviside step function to be zero for time t = 0. This assumption 
excludes the accusal terms in the perturbation theory p6[.In this formalism the correlation and response functions are 
given by path integrals weighted with a density exp{LyFoi: instance if the phonon variable is coupled to an external 
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field h in the entropy functional (2.1), an additional contribution /ik.tdQ-k.-c^ appears in the stochastic functional L. 



Then the phonon response function can be written in the form 

G{k,u;) - (Q-k,-c.gk,c.) = ^J V[S]V[iS]V[Q]V[iQ]V[es]V[ies]V[eL]V[iiL]Q-K-^QK^eML) , (2.8) 

where I?[ ] denotes a suitable integration measure and Z is a normalization factor. The bilinear terms in the Lagrangian 
determine the free response and correlation propagators. Usually, we include the terms proportional to the coupling 
w as well as the term jessL into the interaction part rather then into the free part of Lagrangian. Then, the free 
propagators will be given by the expressions presented in Table |, where c = Cl(^ . Note that in order to obtain the 
free response functions for the spin, spin-energy and lattice- energy fields we have to multiply the response propagator 
by the corresponding damping coefficient . Once the propagators are defined in principle a perturbation expansion 
for any quantity can be performed. Note also that because of the term jesSL, the spin-energy densities es, es and 
the lattice-energy densities cl, sl are dynamically coupled, even if w = in contrast to the static case where they 
are decoupled from each other in the functional H for w — 0. 

C. Decoupling transformations 

We have found it more convenient first to decouple es, es from e^, by the Gaussian transformation: 

■ ' ^ -DLo{k,uj)el^+A{k,uj)Bl^ , 



Csr, 



-DLo{k,-uj)ei 



ClT, 



(2.9) 



Here A{k, oj) 



-DLo{k,uj) 



\DUk,u)\ 



is a coefficient. For the case of the total energy conserved, this 



et,u.Dl^{k,u;)e\_ 



-uTSi^A,-^ - 2fTS\^e\_^ - /A^(fc, w)rs(fc)g^,^5lk,-. + /le^^.^lk,- 
-iwkJ{k,u)Qk,^e1^^_^ - iawkAk-ujQi^.Lje% _^ - ikwK(k,uj)T s{k)e^^^Q^\^^^ 



k,— w 



-iw{aTL{k) - 7)fcek,i^Q-k, 



transformation can be also looked at as a substraction of the slow component associated with thermal conduction mode 
from Ci, so the new variable contains only the fast mode of equalization of temperatures: e^^™ = Cl(^ — ^) 
for Lo = 0. The Lagrangian expressed in new variables takes the form 

[J2 {r^k,^^-k,-c. + dk^Q^^Q-K-^ + rso(fc, w)4c.e^k,-c. + rL(fc)e^^^g^k,-c. 



(2.10) 



where S^k.w = ('5'5')k,w, — {S'^)u,uj and 



Tso{k,Lo)^Tsik)^TL (k) 



-f^\DLo{k,u;) 1^ 



N[k,uj) = l — — ,J(k,uj) = l + 



ClTs 



Cs 



K{k,uj) = N{k,uj) - —(1 — ) . 



As a result of the transformation (2.9) the spin-energy response propagator is dressed to the form 



Dso [k, ^) = (fc, uj) - YsCsCl " ^ 1 



1 ~ iuj 



TsTl (1 - «w) 



(2.11) 
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where uj = u/los and uj = lo/ojl is the frequency reduced with respect to the bare spin-energy = Ts/Cs and 
lattice-energy lul = T^/Cl relaxation rate, respectively. It is easy to see that for very slow lattice-energy relaxation 
^ 1, which takes place for 7^ = 7 ^ or Cl 00, the dressed spin-energy propagator reduces to the bare one 
and the model with two energies reduces to the model without the lattice-energy field, in which only the equation 

SH 

e-s = -^s{k)- h if 

bes 

matters in our considerations. Depending on the ratio of the bare relaxation rate of the energy diffusion mode 
to the bare spin-lattice relaxation frequency 7/C5, the spin system can be then regarded either as thermally isolated 
from the lattice {-^ — ^ 1) or as freely relaxing to an infinite heat bath at each lattice site, or to a bath with an 
infinite thermal conductivity ( y ^1). The same happens for the nonconserved energy case. If 71, ^ 00 the lattice 

energy always takes its equilibrium value. Then w 1, p^p^ — * and again Dsq Dsq. 

The one energy field model we recover also for 7 — > 00 or Cl — > in the conserved total energy case. Then iD <^ 1 

and expanding Eq.( p.l]| ) we get Dso{k,Lu) = Cs^g^m^^ ^1 — i^^) where m{k) = 1 — rs(fclri.(fc) — '^^-^ with 
A_|- = As + Al and C+ = Cs + Cl- being the thermal conductivity and specific-heat of the total S + L system. 
At the same time the dressed damping coefficient Tso {k,uj) rs{k)m so the dressed spin-energy response function 

Tso (k.uj) Dsoik,uj) Cs ^1 — • It can be said that on very short time scales, of an order of Cl/j, the 

lattice comes to equilibrium with the spin system and then for much longer times, of an ord er of the inverse of 
ultrasonic frequency, only the diffusion of the total energy es + bl matters and instead of Eqs.( 2.2c , 2.2(j ) only the 
equation 

e = -A+fc^— + (f> 
oe 

may be used. This limiting case for dynamics was considered by Drossel and Schwabl p^ . The reason for introducing 
the second energy field is now becoming clear. We expect that the model with only one energy density correctly 
describes the dynamics only in the two limiting cases: for pure diffusion of the energy field and for pure relaxation to 
an infinite heat bath. As we want to be able also to describe the dynamics which is intermediate between these two 
cases we have included cl as a separate variable. 

Next, we decouple the sound mode from the spin and energy fiuctuations by the transformation 

Qk,Lj Qk^u - gfcGo(fc,u;)S'k,^ - wkJ{k,u)Coik,uj)e^,_^ - awkCo{k,uj)e^^^ 
+B{k, u;)J\^ + E{k, uj)el^ + F{k, u;)St^ , 

Qk,c. ^ Qk.^ - 25fcGo(fc, -co)rs'\^ - wkK{k, -co)rsGoik, -Lo)el^ - wkiaVL - 7)Go(fc, -c^)e£,^ • (2.12) 
Here the abbreviations 

B{k,uj) = -ATgk\Co\^0k^ , 

E{k, Lu) = -2wkK{k, -u;)Ts\Go\^ek^ - awkA{k, w)Go(fc, uj) , 

F{k,u) = -2wk{arL ~j)\Go\'^ek^ , 
have been introduced. Finally, we apply the transformations: 



+ DLo{k,-Lu) 2awgk'Co{k,-uj)TS\^ + aw^k^K{k,~Lu)Go{k,-uj)Tsii^ , (2.13) 



and 
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+ DsoiK uj)rs [wgk^K{k, Lu)Go{k, cu) - /iV(fc, ou)] S^^ + W{k, uj)S^ 



4,. - eL - 2D5o(fc, [/ - wgk^Jik, ~uj)G^{k, -to)] TS\^ , 
decoupling es and from the spin fluctuations, where 

Dl^ik,Lu) = D^^ik,u;) - (aVL - j)aw^k^Goik,u) , 



(2.14) 



DsQik,Lj) ^ Ds^{k,u;) - w^g^k^Jik,Lj)Kik,u;)Go{k,u;) , 
are "dressed" by transformations ( |2.9| ) and ( |2.12| ) energy-density propagators. Here the coefficients are given by 

Rik, to) = [7/ + (aVL - -i)wgeGo{k, u)] D^k, u^) , 

T{k, uj) = w^ej{k, uj){aTL - l)Go{k, u;)DLo{k, u) , 
U{k,uj)^ATwgk\arL-l)\Gof9k^DLoik,u;)+Aagwk^rrLGoik,-u;)\DLo\^ , 
V{k,u;) = wk \2awkK{k, + (aFi - -f)w^k'\Go\'ee]Go{k, -to)\bLo?T sT l - {oTl - l)E{k, -Lj)DLo{k,u;)^ , 



W{k,u;) = 2wgk^ [aGo{k, -Lu)A{k,Lu) + 2K{k,Lu)rs\Go\^ek^] rDsaik,Lu) 

-4 [/ - wgk'Jik, -uj)Go{k, -Lu)] (^1 - 2:!l^£g^^ rrs\bso\' 



III. ACOUSTIC SELF-ENERGY 



Having apphed the transformations ( p.9| , 2.12-2.14) the acoustic phonon response function, to the leading order in 
coupling constants g and w, can be written as 



(Q-k,-^Qk,c.) = Go{k, to) + Gi{k, c^)E-^(fc, Lo) , 



where 



and 



^rcd 



X{k, uj) = a{aTL - 7)^^ , Y{k, ^ w^TsK{k, Lu)J{k, lu), 



5i(fc,w) = 5 - wfiTsK{k,uj)Dsoik,uj) + awfT^DLoik^uj) 



(3.1) 



{k,Lu) = fc2 {x(A:,L^)^Lo(fc,t^) +l'(fc,c^)^5o(fc,t^) + 2gi{k,Lu)g2ik,u;){TS^-k,-u.SlJ''^"} , (3.2) 



g2(fc,w) = 5 - wf2Jik,uj)TsDsaik,uj) , 

with / = 7/ + (arL_7)w.gfc^Go(A:,t^), /i = fN{k,uj)-wgPK{k,Lu)Goik,Lu) and /2 = f -wgk^J{k,uj)Goik,iu). The 
crucial point in Eq. (pj) is that the four-spin response function is calculated with the Lagrangian which depends solely 
on the spin variables. As the result of coupling of the spin varia ble to the energy densities and phonons, the effective 
spin Lagrangian Lcs obtained by transformations ( p. 9 , 2.12| -2.14) contains phonon- as well as energy-density-mediated 
four-spin non-local interactions u{k, uj)S^ -k,-uiS^ ^ with 
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u(fc, Lo)^u- 2g^k''Go{k, to) - 2awgk'' fVLGoik, u;)^Lo(fc, uj) - 2hf2TsDso{k, u) . (3.3) 
A general expression for the interacting phonon response function is given by 

G-^{k,uj) = Go^{k,oj)-Y,{k,oj), (3.4) 

where the phonon self-energy E(fc,Ci;) is the irreducible, with respect to phonon lines, part of Y7'"^{k^uj). From Eqs. 



(3.1) and (|3.4|) the irreducible and reducible parts are related through 



^(^'")=l + Go^fc,i)E-d(fc,^)- (3.5) 

Next we may eliminate the dangerous resonances [|2|,pj by replacing Go{k,uj) with G'o(fc,0) in S(fc,w) i.e. we set the 
strongly irrelevant parameters - the coefhcients at wfc and oj^ - equal to zero. In this paper we will not discuss the 



role of the denominator in Eq. (3.5) assuming the elastic couplings g and w to be very small. This assumption gives 



I](fc,w) {k,Lu). It is consistent with neglecting the higher order (in g and w) terms in Yi {k,Lu). Then also 

Z)io(fc,^) — DLQ{k,uj) and Dso{k,uj) = Dgo(k,uj) together with small g and w let us also neglect the macroscopic 
instability P,|l8t , which is believed to take place in compressible spin systems with positive specific-heat exponent near 
the transition temperature. The weak first-order transition in such systems is a result of a non-analytical character 
of the coupling constant u{k,uj) at A; = 0. For small g and w the first-order regime can be probably observed 



only extremely close to Tc- We could in principle put other irrelevant parameter in Eq. (5.5) equal to zero as the 
coeflScient at lo in the energy-density propagators for non-conserved systems but then we would obtain only the 
asymptotic behavior of the sound attenuation coefficient. As we are also interested in the nonasymptotic behavior 
associated with the coupling to the energy-density we will proceed in another way, preserving the nonasymptotic 
effects. First, note, that because for small g and w 

u{k, u;) = u- 2fN{k, Lu)TsDso{k, lu), 

thus (rS'^_k,-(^5'^ uj)^"" reducible with respect to (new) spin-energy propagators and so is S . We assume that the 
coupling constant / need not to be very small but the positivity of the four-spin coupling must be guaranteed [||. 
The perturbation expansion with respect to u{k,uj) let us express (r5^_k,-t»j'5'k i^)^"", the reducible with respect to 
Ds four-spin response function, in terms of the irreducible one 

(rs^-k,-.^^ J"- = {rs^-K-^slJ^'^ ^^^^^ 

1 - 2/2rsA^(fc, uj)Dsoik, c^)(r52_k,-^^^ J'" 



Eq. (3.6) is equivalent to 



1 + 2prsNik, u;)Dso{k, u;){rS^^^^^^Sl Ji^^« 



It is easy to note that the coefficient preceding (TS'^ -u,-u)S'^ ^)^^" in the denominator of Eq. (3.7) can be written 



as 2pDs{k,uj) with Ds{k,uj) = Gs {I — iujb{q,uj))~^ where b{q,uj) = ^^zj^- -'-'^ a; = fimit it is equal to 
V = 2pGs- As (TS'^-]s_-^Sl^^Y" is a vertex function, containing no resonances, we may again replace strongly 



irrelevant parameters in Eq. (3.7) like the coefficient at uj in Dso, in the case of nonconserved systems, by zero. We 
then obtain 

= , (3.8) 

where La is the action of the model A of Halperin, Hohenberg and Ma ||l^, with ua — u — v. In the case of a 
system with the conserved total energy we cannot proceed in this way as Dso contains nonstatic relevant terms 
and (rS'2_k,-i^S'^ jj)"''^ should be calculated using energy generated dynamic interactions in L^s [1). The asymptotic 
behavior of the ultrasonic attenuation coefficient in a model with conserved energy field was investigated by Drossel 
and Schwabl . They showed that the difi'erence of the asymptotic scaling functions (compared to the model A) 
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originates mainly from different exponent z for model C. In the present paper we will restrict our discussion to the 
nonconserved systems an d a detailed analysis of the co nse rved systems will be presented elsewhere. 



Now we return to Eq. (3.8) which inserted into Eq. (|3.6|) yields 



^2 \Le" (r5'2_k,-(^'S'k,^) 



(r52^k,_^5^ )^e« = ^ . (3.9) 



Next from Eqs. (^,(|3j) and (|3j) we obtain 

A ~ 25i(fc,w)52(fc,a;)(r52_k-.^S'2 )iA ] 

1 l + 2/2(C5-i?s)(r52_k,-.^U^M 

It is not difficult to note that 

gi{k,oj) = g2{k,oj) = g{k,uj) = g- wfP{k,uj)Ds{k,oj), 

with 

P(k,Lu) = ._ . 



Now let us split ADlq + BDso into the following terms 



■u?k^P'^Ds{k,ijj) 



and reduce the second term to the common denominator with the third term in Eq. (3.10). This results in the 
following expression for S 

I](fc, uo) _ a^w^CLm , [Iw^^P^CsDsik, uj) - 4wfgPDs{k, uj) + 2g^] {TS^^ ^^.^^ Sl J^^ + P^ 

. (•J-llj 



l + 2p{Cs-Ds){rS^-i,,-^S, 



1 \L. 



Remembering that Cg — Dgik^io) = —iCijb(k,Lo)Ds(k,io) and ignoring the first noncritical term in Eq. (3.11) the 
acoustic self-energy can be written as 

S(fc,c.) ^ 2{g\K^)~i(l:h{k,uj)g^\l,vJ^^-^^^SlJ^-+w^P^Cs ^2) 
l-z(i6(A:,c^)(l + w(r52_k,-^S'2,^)^A) 

with g{k, uj) = g — wfC'sP{k, uj). 

IV. DISCUSSION 
A. General expressions 



The critical contribution to the coefficient of attenuation is determined by the imaginary part of S and Eq. (3.12) 
implies 

2c-^ |1 — ia;o(l + tiw)p 

where * = (rS^_k,-c^S'^.^)-^* and Wi = 2vQRe[bg^{k,uj)], 

W2 = 2Rcg^{k,Lu) + 2Q^\b\^g^ +Qlm[vbw^CsP^ + 2bg^ - 2bg^{k,Lu)], 

W3 = 2CjRe[bg^{k,Lu) ~ bg^ + bw^fClP^] + 2ling^{k,Lu). 
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In the limit — > cxd these coefficients take the corresponding vahies from the model without lattice-energy field Q 
i.e. Wi = 2vojg'^, W2 = 2(g^ + w^g^) and W3 = —AwfCs(^g with g = g — wfCs- The four-spin response function 
usually is evaluated in the limit fc = since in the ultrasonic experiments the wavelength is much longer than the 
correlation length (fc^ ^1) while the ultrasonic frequency can be comparable to the characteristic frequency of the 
spin fluctuations. Then the singular part of the four-spin response function \l/ obeys the scaling relation 

* = t-'^Hy) 

with y = ujt~'"^ /T as the reduced frequency and t proportional to the reduced temperature. Here a, v and z are usual 
critical exponents. The scaling function $ is known to the leading order in e = 4 — d : 



$(y) = Q-'^'^ I - + - 
la y 



I (1 - iy/2) arctan(y/2) - - In (l + {yj^f) 



where 6 = [1 (y/2)2]-4 and = (STr^)"^ With the aid of Eq. (O) the attenuation coefficient can be written as 



-^Wir^'^F^iy; w, t) + -#2i"^"+"'^'F2(y; t) + W^r'^F^iy- uj, t) 
1 5 1 



(4.2) 



where Wi = WiK^/Q, W2 = W2K4, W3 = W^K^/uj and Fi = |$/Afp, F2 = Im{^)y-^\M\-^ , F3 = Re($)|Af|-2 
with the denominator of Eq.(3.12) denoted by 



M{t, uj) = 1- iub{l + vt-°'^). 



B. Regime I 

Eq. ( |4.2| ) shows many different regimes. First let us consider the regime I: m > (cj < ) for which 

P{uj) ^ 1. It seems that the regime I may be adequate for some metals, where should be understood as the spin 
energy of the conduction electrons and also infinite specific-heat associated with phonons can be assumed. 



1. Low- frequency region 



In the low-frequency region, ~ ^ ^ 1, the denominator in Eq. (4J) does not play any role and the second 
term in the numerator with strong singularity dominates 



a{uj,t) 



9^ 2,- 



Ini$(y) 



V 



(4.3) 



This type of singularity in the ultrasonic attenuation, which we shall also call Murata-Iro-Schwabl behavior, has been 
so far obtained by neglecting the energy density fields and it is believed to take place in magnetic metals 
Note, however that the coefficient W2 in this limit is equal 2{g — wfCsY to be compared with 2g^ in the model 
without coupling to the energy fields Q| . 

When — ^ 1 i.e. for M {t, w) ~ 1 while t is not extremely small there may be a competition between the first and 



the second term in Eq. (12). For t > t, 



a{u!, t) 



-g'vCs 



IS 

Csr 



u;2i-2"|$(y)|2. 



the weak-singularity term dominates and 

(4.4) 



Such behavior was first obtained by Kawasaki [|7[|. He postulated phonon-spin-energy coupling in order to explain 
extremely small sound attenuation exponents observed in magnetic insulators. Note however the difference between 
the corresponding scaling functions. What also arises from our analysis, is that such behavior is generated even if 
there is no direct coup ling of the sound mode to the energy-density fields, and that it cannot be truly asymptotic as 
both terms ( |4.3| ) and (4.4) are proportional (in the w ^ limit) to the square of the same coupling constant g so 
even if 5 = the strong singularity described by Eq. (4.3) is still present and become dominant for sufficiently small 
t. The dominance of a given type of behavior does not depend on the relative strength of the coupling constants g 
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and w but rather on the ratio of the bare relaxation times and also on the ratio of v to its fixed point value 
V* ~ {a/ u)K^^ . Slow spin-lattice relaxation (and/or fast decay of spin fluctuations) together with strong coupling v 
favors the weak-singularity type behavior. 

It should be noted that the amplitudes of the Murata-Iro-Schwabl as well as the Kawasaki term are proportional 
to the same coupling constant only in the limit w = (D = 0. For finite w or w both amplitudes begin to differ 
from each other and become frequency dependent as can be seen from the expressions for the coefficients VFiand W2- 
For example Wi is not longer equal to W2 — 2Cj'^\b\'^ . To be exact this frequency-dependence (except 2Cj^\b\'^g'^ 
term in W2) of the coefficients Wi results from the inclusion of the lattice-energy field to the model in which spins, 
phonons and spin-energy modes interact. The nonuniversal frequency-dependent coefficients may be useful in making 
comparisons with experimental results. 



2. High-frequency region 



As the sound frequency increases, uj/m approaches unity and the denominator in Eq. (4.1) as well as the second 
term in the expression for W2 become important. With a further increase of the sound frequency, two scenarios, 
can be realized. In the first one, W2 <IC Wi i.e. the sound frequency is not high enough and/or g is much smaller 
then g. Then, the situation resembles the high frequency regime for Kawasaki classical relaxation function The 
sound attenuation coefficient saturates and Ci{ui, t) cx const. The second scenario happens when W2 in which the term 
25^0;^ dominates, is much greater then Wi and then: 

(i) for vt^°'<^ ^ 1 i.e. with noticeable singularity in the specific-heat a new type of behavior can be observed 

a(c.,0oc5W(— )^^^^^^, (4.5) 



which can be obtained from the one described be Eq. (4.3) by a simple replacement a — > —a, $ ^ $ ^ and g 



giving a bit weaker than Murata-Iro-Schwabl singularity but much stronger than the anomaly in the energy-relaxation- 



dominated region described by Eq.(4.4). The difference in exponents zv -\- a and zv — a \s very small and it may be 
difficult to distinguish both regimes in experiment. This new type of behavior of ultrasonic attenuation in magnets 
is very interesting. In this regime the frequency of the sound mode is much higher than the (modified b y th e lattice) 



spin-energy relaxation frequency so contrary to the low-frequency region described by Eqs. (4.3) and {iA) now the 
"lattice" and the spin systems no longer can be treated as remaining at a local equilibrium. The local temperature of 
the spin system is not able to follow alternate hot and cold temperature variations produced by the ultrasonic wave. 
One can say that we have here some kind of "adiabatic" sound propagation although the "lattice" (whatever it is in 
regime I), in virtue of the assumption <C to < 1 defining the region I, can be treated as at equilibrium with the 
ultrasonic wave. 



The behavior described by Eq. (4.3) can be regarded as dominant for t he w hole asymptotic region t, w — > except 
only for a very narrow area vt^°'u}/m ^ 1 where the "adiabatic" behavior (^^) is dominant. The latter area, however, 
can be experimentally inaccessible as for sufficiently low sound-frequency the reduced temperature needed to make 
vt~°'Lb/m big enough may be too small, as a consequence of smallness of the exponent a, to be approached in real 
experiment. For this reason we shall sometimes call the Murata-Iro-Schwabl behavior also the asymptotic behavior 
keeping in mind that it concerns only the experimentally accessible part of the asymptotic limit. On the other hand, 
in many ultrasonic experiments the measured frequency range is limited only to one or two decades so for the magnets 
with slow spin-lattice relaxation this asymptotic region can not be approached at all. Because the "adiabatic" behavior 
comes from the 2uP'\h\'^g'^ term in the coefficient W2 we can also say that as the reduced frequency Cj increases, the 
Murata-Iro-Schwabl behavior turns into the "adiabatic" o ne. 

(ii) For vt~°'^ <C 1 we recover the asymptotic behavior (O) but with cf^ replaced by g^ . What is rather unexpected 



we observe the asymptotic behavior in quite a nonasymptotic region because here (D ^ 1 and t is also large. Moreover, 
the proportionality coefficient g^ is exactly equal to that from the theory which neglects all effects of the energy- 
density fields. However this region is not fully critical because as a result of smallness of vt~°'^ also the critical 
contribution to the specific-heat is very small so the specific-heat does not show any noticeable divergence in this 
region. The last statement is easily seen in the static case, where the temperature dependent specific-heat is given by 
Cs{t) = Cs(l -I- f ('S'^'S'^k)ff)' remembering that the four-spin response function 2(rS'2_k ^.^S*^ transforms into 
{S^S'^^^)h for Lu = 0. Here the static average {...)h was calculated with the weight exp{—H). 

In the above we have assumed that despite of a) S> 1 the frequency to is still small. This assumption not always 
seems to be correct and if so the reduced frequency iD may become comparable with 1. Then the coefficients Wi and 
b{uj) began to evolve, however the "adiabatic" limit (J) — > 00) does not change as the term 2(1)^151^5^ iii ^2 is reduced 



with the dominant term in the denominator of Eq.(4.1) 
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3. Comparison of the "adiabatic" limit with critical ultrasonic attenuation 
in binary mixtures and liquid He 



The analytic relation (4.5) resembles the corresponding asymptotic formulas for sound attenuati on in binary mix- 

. Looking at Eq.(4.5) it is easy to 



tures |19| , |2Cl| , pi| , [L2[ , where the sound attenuation exponent is also equal to zi/ 
see that it can also be transformed into the form 



a{uj,t) oc — ,g^wlni(l + t;(rS'2_k,-ij'S'k.< 



(4.6) 



in this region, which is also identical in form to the expression used by Ferrel and Bhattacharjee p^ ] for the ultrasonic 
attenuation near the A-transition in liquid helium (although the specific-heat exponent is very close to zero for such 
systems) : 



a{uj,t) oc — LjImCi?B(i, t^) 



(4.7) 



where CpBit:^) is a phenomenological frequency-dependent specific-heat. Comparing Eqs.(4.6) and (4.7) we get the 
correspondence between these two equations by interpreting 



(4.8) 



as a frequency-dependent specific-heat in the Ferrel-Bhattacharjee sense, where as usual we put fc = 0. A very similar 
identification, in the case of ^He, was made by Pankert and Dohm who gave the statistical meaning to CpBiti^^)- 
Note however that in spite of the fact that for w ^ Eq.( 4.1 
and g) into the exact relation known from the statics 



transforms (up to bilinear terms in small couplings w 



Cs{t,k) ^Cs + {f- wgcl)'Cl{SlS\)H, (4.9) 
where Cs{t,k) = {e^e^]s}H is the static temperature and wave vector dependent specific-heat, with Cg^ — C 



uq and Cq 



-1 
s 



Cl, the interpretation of Eq.(46) in terms of true frequency-dependent specific- heat is 



inaccurate as in dynamics the latter is given by the product of the spin-energy propagator Ds{k, lo) and the frequency- 
dependent damping coefficient Ts{k,u!), which can be found as the coefficient staying in front of the egeg term in the 
Lagrangian after transformation ( |2.14| ). From the transformations (2. £ ,2.13-2.14) we find 



Ds{k, uj) = Dso{k, uj) + 2^5o(fc, iu)h{uj)f2{u;){rS^^ 



\ioft 



(4.10) 



and Ts{k, uj) ~ rso(fc, uj) if we drop the terms of o rder of w'^. The determined in this way frequency-dependent specific- 
heat, although giving the same static limit (O) for w — > differs significantly from Cg{t,uj) for high frequencies 

{lu ^ oo) where Dsoik, uj) ^ contrary to Cs which is a constant. Also the four-spin response functions in Eqs.(^?^) 
and (4.10) show different behavior as the coupling constant u is frequency independent in the idealized Lagrangian 
La in contrast to u{k,uj) in Loff (Eq.(3.3)). The advantage of the function C^(i,a;) is that it is much simpler in 
calculations, as La does not contain the details of the dynamics. It is worth noting that C'^{t,uj) is defined with the 
help of the constant v = 2f^Cs which is assumed to be finite even if the coupling constants g and w are very small 



i.e. also in the so-called "we ak- coupling" limit [pUlS]]. 

The "adiabatic" formula (4.6) shows a close analogy to the critical sound attenuation in two quite different systems 
belonging to different universality classes and in the case of *He even the order parameter dimensionality is different, 
n = 2 and a ~ 0. It may suggest that as concerns the adiabatic sound propagation this kind of singularity may 
be quite common in the critical systems where the order parameter is coupled to the energy and sound mode by 
two different coupling constants. The exception is the gas-liquid system but there, as was noted by Pankert and 
Dohm p^ ] , the order-parameter itself is proportional to the energy and there is only one static coupling between the 
order-parameter and the sound variable. 



C. Regime II 

Now let us consider a regime which is more appropriate for description of the common lattice consisting mostly of 
phonons. The new regime, opposite to regime I, is defined by the inequality: m <^ lj {uj ^ — — -^^^-^—)- Now the 
fimction b{uj) -Naj(i-m) ^ ^^^^ assume further that m c± <C 1 and m ^ ui^ then b{uj) ^ 1 + ^, 
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P{lo) ~ 1 — a + ia^ and —iCohiuj) ~ ^(1 — iuj). Inserting these expressions into Eq.(3.12) and neglecting the terms 
of order of ™, or higher, we get 



vCsnk.^) _ 2g^CLo{to)C§{t,uj) + 2g^CsAC§it,uj) + [w^^ - a)^Cs - 2g^]CsCLo{^) 



fc2 



(4.11) 



where Clo(^) 



is the free, frequency-dependent, lattice specific-heat, g = g — wf{l — a)Cs and ACg 



C-^{t,uj) — Cs is the fluctuation only contribution to C-^{t,uj). 



1. Low-frequency and high-frequency regions 



It is seen from Eq.( 4.1l| ) that as long as | Choi'-o) |^| Cg(t,Lu) | i.e. when u; is not too high the denominator is only 



slightly singular and we have a kind of competition, with frequency-dependent weights, between two kinds of low- 
frequency contributions, the Murata-Iro-Schwa bl a nd Ka wasa ki term. Note only, the different amplitude, proportional 



to g^ in the limit d) = 0, instead of 5^ in Eqs. (4.3) and (4.4). In this low-frequency region again the spin and lattice 

Cs, which is equivalent to the 



Cl 



systems stay at equilibrium with each other. For very high tu: \ Clo{'-^) 
condition u) » 1 the specific-heat C^{t,uj) begins to outweigh in the denominator of Eq.( 4.11 ) and the "adiabatic" 
term proportional to 



— const 



g'Cs 



1 



dominates the acoustic self-energy, because the spin and lattice subsystems are no longer in equilibrium. 



2. Relevant frequency-dependent specific-heats 
Introducing C^{t,Lu) = and C^{t,Lu) = Clo(w) + C^{t,i^), Eq.(pl|) can also be written as 



vCsm,^) _^^2u ^r^A., , , 2g^CsAC§it,co) + [wMl-arCs-2g']CsCLo{^) 

p ~ 25 {t,u;)C_{t,u) + . (4.12) 

C^{t, oj) is a kind of specific heat which by analogy with the static, lu = 0, case can be given the following interpretation. 
The energy contribution to the functional H, describing the statics, can be rewritten as 

22 22 
Cs Cl C+ C_ 

where 1/C_ = I/C5 + 1/Cl and ca = C_(^ ~ cT-'' ^^^^^^ lattice and the spin systems are at equilibrium with each 
other, then the local fluctuations of temperature 6Ts = ^ and ST^ — ^ are equal, and consequently A = eA = 0. 
Thus CA can be interpreted as excess energy in one system (and deficiency in the other) in comparison with the case 
when they were at equilibrium with each other {STs = ST^), so if we shift the energy SEs from the spin system to 
the lattice system {SEl = —SEs) it results in the shifts of temperatures STs = and STl = = where 
V is the total volume, in both subsystems. The total temperature difference 



SEs SEs SEs 

Ts~n^STs-Sn^^ + ^ = ^ (4.13) 

is proportional to the energy transported from one system to the other and C-V — Xs-Tl understood as the 

net heat capacity for the process of differentiation (or equalization) of temperatures between both subsystems, whereas 
the total specific heat C+ ch aracterizes the thermal properties of the system as the whole (after equilibrization). 



Coming back to Eq.( 4.12 ) note that it can also be rewritten as 



vCs^jk^co) 2 -J, 2{g^~g^)CsACi{t,co)+wM^-a)^ClCLo{u^) 

2g (t,^)C_(t,^) + (4.14) 
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where C^{t,uj) = '^^Lot'^'i+Cs+^Ht^) looked upon again as the net specific- heat for the process of equilib- 

rization for two sub-systems with the only difference that this time the background term Cs has been "moved" from 
the spin to the lattice sub-system. Now the lattice degrees of freedom plus short-wavelength spin fluctuations consist 
toget her t he one (noncritical) sub-system and the long- wavelength (critical) spin fluctuations the other. In Eqs. ( 4.12| ) 



and ( 4.14 ) the second terms vanish for g = and w = Q, respectively. We shall call them the mixing terms because in 
the case when one of the coupling constants g or wf{l — a)Cs dominates the acoustic self-energy is composed of the 
dominant term proportional to the corresponding specific- heat C^{t,u) or C^{t,uj) and a small correction resulting 
from a slightly different ways the constants g and w couple the sound to the spin system | p5[ | . The former one couples 
the sound mode only to the part of the spin system, the long wavelength critical spin fluctuations, in contrast to the 
latter coupling which couples effectively the whole spin system (including the noncritical background) to the sound. 
Looking at C^{t, cu) it is easy to see that for Cl ^ Cs it contains both low- frequency terms: the Murata-Iro-Schwabl 
and Kawasaki terms, as well as the "adiabatic" limit. Namely, the first term in 

I AC§{t,co) 1^ lmCLo{io)+ I CLojco) + Cs \' lmACi{t,Lo) 
I Clo(w) + Cg[t,Lu) \^ 

can be identified as the Kawasaki term for a) ^ 1 {CLoiuj) ^ Cg(i,cj)) and the second as the Murata-Iro-Schwabl 

one. The crossover temperature Across — (ilcT^J^^cIr) depends mostly on the ratio of bare relaxation frequencies 

ws/r. In the opposite case a) » 1, CLoi^) ^ Cg{t,u)), and ImC'^(i,w) ~ '"^ica^^^)]*'""' i-^- "adiabatic" limit is 
recovered. The situation described above may be also superimposed by the effects associated with the crossover in 
the specific-heat Cg{t,uj), because the expression for Cg{t,ui) — Cs(l -I- wi~"<I>) may not display the required 
singularity in the physical region of interest (if v is very small). Then, for example, the measured effective sound 
attenuation exponent in the Kawasaki region will be much smaller then 2a, sometimes close to zero. 

In closing we shoul d men tion that we have not discussed yet the third term in the numerator of Eq. ([1.2|) as well 



as its analogue in Eq.(4.11). It gives the weakest singularity oc which is usually smaller then the first two terms. 
However, as it is proportional, in the w = limit, to wfCsg i.e. to the first power of g, whereas Wi and 11^2 are 
proportional to the second power, then the i^" term may dominate the ultrasonic attenuation in the degenerated 
case of 5 — > 0. 



V. SUMMARY 



We have performed a detailed analysis of the critical sound attenuation in magnets. The relaxation of spin-energy 
to the lattice has been fully taken into account in the prediction of the temperature and frequency dependencies of the 
acoustic self-energy. The important point here is that the acoustic self-energy is reducible with respect to the energy 
propagators. We were able to express Y,{k,uj) in terms of the four-spin response function {TS'^ -■^^-coS'^ ^)^^ of the 
idealized, phonon- and energy-free model A. The last quantity can be relatively easy calculated by the renormalization 
group method. In the low-frequency region, two kinds of critical behavior described by Eqs. (^^) and (4.4) compete 
with each other and the relative weights are determined by the coupling constants g and w as well as the frequencies 
w, ws, and the other parameters. The high-frequency or "adiabatic" limit is determined completely only by one 
coupling constant g - describing the interaction of sound mode with two spin fluctuations. The dynamic scaling 
functions for these three regimes are calculated within the one-loop approximation. 

Another point of interest is the possibility of expressing the acoustic self-energy by a suitable frequency-dependent 
specific-heat. We have shown that in the case of one coupling g or wfCs{l — a) being much stronger then the other, 
the acoustic self-energy can be very well approximated by the specific-heat C^(i, uj) or C^{t, lu), the net specific-heats 
for the process of equilibrization between the lattice and spin sub-systems or in the latter case between the noncritical 
background composed of the lattice degrees of freedom plus the noncritical short-wavelength spin fluctuations and the 
critical subsystem composed of the critical long-wavelength spin fluctuations. The difference in deflnitions of C^{t, uj) 
and C^{t, uj) comes from a slightly different ways the couplings g and wfCs{l — a) couple the acoustic phonon to the 
fluctuations of the order parameter. The virtue of Eqs. (4. 12) and (4.14) is that (depending on the relative strength 
of the coupling constants) in some regions of parameters the acoustic self-energy can be expressed entirely by the 
frequency dependent speciflc-heat C^{t,ijj) (or C^(t, a;)), governing the equalization of temperatures between the 
lattice and spin systems. C^{t,uj) shows its relevance in the ultrasonic propagation in magnetic systems due to the 
fact that depending on the frequency, the reduced temperature and the ratio ws/T it can reveal the most important 
types of singularities in the critical sound attenuation. 
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A similar analy sis is also possible for the magnets with the order parameter dimensionality greater than one. For 
example the Eq. ( 3. IS ) is still valid there. However, the specific- heat exponent a is negative in X-Y and Heisenberg 
systems, therefore, the scaling behavior of the Kawasaki as well as the "adiabatic" limit will change. 
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The mixing term in Eq.(4.12) participates also in forming the correct "adiabatic" limit which always is proportional to g 



TABLE I. Propagators of the model 



Propagator 



Analytic expression 



Free Lagrangian average 



Spin response- 
Spin correlation- 
Phonon response- 
Phonon correlation- 
Spin- energy response- 

Spin-energy correlation- 
Lattice-energy response- 

Lattice-energy correlation- 



Kso(k,iL!) = 
Go{k,uj) = ^^j^ 
KQo{k,Lo) = 
DsQik,Lj) 



icj + r(r + fc^) 

1 



1 



Keso{k,ijj) 

DLo{k,Lu) = 



- ,2rs(fc) 
' ^^+rs(fc)^ 
1 

-ic^+ri,(fc)C7 



(>S'k,t^S'_k,-ij)o 

{Sk,UjS-k,-Lj)0 

(Qk,t^Q-k,-cj)o 

(ek,"e-k,-c.))o 
(ek,c.e?k,-c.))o 



.))o 
.))o 
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